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Neutrons and antiprotons in ultrahigh energy cosmic rays 
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Abstract. The neutron fraction in the very high energy cosmic rays near the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff energy is analyzed by taking into account the time dilation effect of the neutron decays 
and also the pion photoproduction behaviors above the GZK cutoff. We predict a non-trivial neutron 
fraction above the GZK cutoff and a negligibly small neutron fraction below. However, there should be a 
large antiproton fraction in the high energy cosmic rays below the GZK cutoff in several existing models 
for the observed cosmic-ray events above and near the GZK cutoff. Such a large antiproton fraction can 
manifest itself by the muon charge ratio n + / fjT in the collisions of the primary nucleon cosmic rays with 
the atmosphere, if there is no neutron contribution. We suggest to use the muon charge ratio as one of the 
information to detect the composition of the primary cosmic rays near or below the GZK cutoff. 



The observation of ultrahigh energy cosmic rays above 
and near the Greisen-Zatsepin-Kuzmin (GZK) cutoff en- 
ergy presents an outstanding puzzle in astrophysics and 
cosmology 1 . It has long been anticipated that the high- 
est energy cosmic rays would be protons from outside the 
galaxy, and there is an upper limit of the highest energy 
in the observed proton spectrum, commonly referred to as 
the GZK cutoff [2], as the protons travelling from inter- 
galactic distances should experience energy losses owing to 
pion productions by the photons in the cosmic background 
radiation. Then it is suggested [3] that the particles with 
energy above the GZK cutoff must be coming from within 
the local "supercluster" of galaxies of which we are a part. 
Thus, the "GZK cutoff" is not a true cutoff, but a sup- 
pression of the ultrahigh energy cosmic ray flux owing to 
a limitation of the propagation distance, which we refer to 
as the GZK zone. Although there have been some novel 
explanations for the observed cosmic ray events above the 
GZK cutoff, it is natural to expect that these ultrahigh 
energy cosmic rays come from sources within the GZK 
zone, i.e., not far from us in more than tens of Mpc. The 
"Z-bursts" [lEE] can explain the highest energy cosmic 
ray events, with sources within the GZK zone [SJ- The 
reason is that the "Z-bursts" are from the Z-boson anni- 
hilations of the ultrahigh energy neutrino (antineutrino) 
cosmic rays with the relic neutrinos (antineutrinos) in the 
cosmic background radiation, and these annihilations can 
happen anywhere in the universe. The "Z-bursts" can pro- 
duce nucleon cosmic rays with ultrahigh energy within the 
GZK zone, as the energy of the produced Z-bosons is high 
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enough by the collision of ultrahigh energy neutrino beams 
with the relic neutrinos of non-zero mass. Therefore the 
nucleons with ultrahigh energy from "Z-bursts" can reach 
us without the GZK constraint, and this can explain the 
observed cosmic ray events with energy above and near 
the GZK cutoff. Several existing models of similar kinds, 
such as using relic metastable superheavy particles fJJ and 
topological defects |S] as origins of ultrahigh energy par- 
ticles, can also explain the ultrahigh energy cosmic ray 
events above and near the GZK cutoff. 

In fact, these models can produce a number of pho- 
tons and nucleons to reach the earth, for example, the 
photon/nucleon ratio is of about 10 on average in the Z- 
burst model. A possible method to test these models is 
by the large j/p ratio for the cosmic rays near the upper 
end of the spectrum. This can be achieved by the sepa- 
ration of primary photon and proton cosmic ray events 
with different characteristic air shower profiles. However, 
for the nucleon cosmic rays from these models, there could 
be a sizable number of neutrons and antinucleons in com- 
parison with that of the protons. For example, there are 
roughly equally numbers of protons and neutrons at the 
source point in the Z-burst model, with also a symmetric 
part of antinucleons. Therefore it is necessary to estimate 
the composition of the nucleon cosmic rays. We will study 
the role played by the neutrons in the ultrahigh energy 
cosmic rays near the GZK cutoff. We will take into ac- 
count the the time dilation effect of the neutron decays 
and the pion photoproduction behaviors of the nucleons 
above the GZK cutoff. It will be shown that there should 
be a non-trivial fraction of neutrons above the GZK cut- 
off, whereas the neutron fraction is negligibly small below 
the GZK cutoff. Since there should be equally numbers 
of protons and antiprotons in the cosmic rays from the 
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Z-bursts, and these protons and antiprotons will produce 
different muon charge ratio fi + / fi~ in the air showers, we 
thus suggest to measure the fi + / fx~ ratio of the air showers 
by the primary cosmic rays near or below the GZK cutoff, 
as one of the information to detect the composition of the 
primary cosmic rays with high energy. 

The 2.73 K cosmic microwave background (CMB) of 
the photons satisfies the Planck's ideal black-body radia- 
tion formula, with the number density n 1 = 167r£(3) (kT /he) 
413 photons per cm 3 and the mean energy per photon 
e 7 = 7r 4 fcT/30C(3) = 6.35 x 1CT 4 eV, where C(3) = 1.20 
is the Riemann Zeta function. When the nucleon with 4- 
momentum p = (E, p) interacts with the photon with 
4-momentum k — (e, k), and composes into a system with 
the center of mass energy squared S, we have 



n + Jcmb —> Ti" + n. These four processes have quite dif- 
ferent cross sections, as can be understood from the low 
energy theorem [S] as well as from the chiral quark model 
[1T)j . Roughly speaking, the cross section of the pion pho- 
toproduction N(j, ir)N is proportional to A 2 with the four 
physical amplitudes of A expressed by [TU] 



E=(S 



/2c 1- a/1 ~m 2 N /E 2 cos 



(1) 



where is the angle between p and k. cannot be zero 
since a nucleon cannot catch up a photon moving in the 
same direction, and the energy of the nucleon E must 
be very large near the pion photoproduction process N + 
Jcmb —* 7T + N, therefore we have, 



E sa (S 
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The threshold energy for pion production N + Jcm b 
7T + N is 

E w (2m N m 7r + ml) /4e = 1.10 x 10 20 eV, (3) 
and the threshold energy for producing the A resonance 

+ JCMB -> A -> 7T + N is 



E 



I 2 2 

\m A — m 



N 



) /4e = 2.52 x 10 20 eV 



(4) 



The neutron has a mean life time r n = 887 s in its rest ref- 
erence frame. Due to the time dilation effect of the special 
relativity, the life time of a moving particle is dilated by a 
factor j n = E n /m n . Therefore the mean free path that a 
moving neutron can travel before its beta decay into the 
proton should be, for a neutron at the pion production 
threshold, 



l n ~ c7„t„ = 3.12 x 10 24 cm = 1.01 Mpc, 
and for a neutron at the A resonance threshold, 
l n ~ c lnT„ = 7.11 x 10 24 cm = 2.30 Mpc. 



(5) 
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Therefore the neutron fraction around the GZK cutoff is 
expected to be negligibly small for nucleons coming from 
a source with distance of considerably more than a few 
Mpc, if we only consider the time dilation effect of the 
neutrons. 

However, the situation is different if we take into ac- 
count the detailed features for the pion photoproduction 
N + Jcmb - * I" + N, which, when distinguishing between 
proton and neutron, should be written as p + Jcmb — * 



A(pj -» 7T+n) =V2(A~ + A ) 
A(nj -► tt~p) =V2(A~- A ) 
A(pj->n°p) = A + +A°, 
A{nj -> 7r°n) = A + - A , 



(7) 
(8) 
(9) 
(10) 



where the isospin amplitudes can be expanded in terms 
of the ratio between pion mass and nucleon mass r\ = 



A- = l + 0(rj 2 ), A + = A = -r)/2 + 0(r? 2 
Thus we get 

a(nj —> tt~p) (I + 77/2) 2 



a(pj 



-n) (I-77/2) 2 



1.34, 



(11) 
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which is in excellent agreement with the experimental data 
HB, and 
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0.01, 



(13) 
(14) 



which means that the neutral pion production processes, 
P + Jcmb — ► tt° + p and n + Jcmb — * n° + n, can 
be neglected. Adopting an average cross section <j(pj — > 
) = 200 fib 0E] above the pion photoproduction 
threshold, we have the mean free path of interaction for 
the proton 



1 



1.21 x 10 25 cm = 3.92 Mpc, 



n 1 a(pj — > 7r+n) 



and that for the neutron 



(15) 



9.04 x 10 24 cm = 2.93 Mpc. 



y a{nj 



-P) 



n, n + jcmb -> tt + P, P + Jcmb 



p, and 



(16) 

It is interesting to notice that the protons and neu- 
trons change into each other via charged pion production 
by the relic photons in the travel until the nucleon ener- 
gies degraded to below the GZK cutoff. There is always a 
certain amount of neutron fraction for nucleons above the 
pion photoproduction threshold, since the protons can al- 
ways change into neutrons via the charged pion photopro- 
duction. Though the neutrons change more fast into pro- 
tons via both beta decay or charged pion photoproduction 
(with effective mean free path A° ff = l n \ n l (ln+^n)), these 
produced protons continue to change into neutrons if their 
energy is still above the pion production threshold. As a 
consequence, there is always a non-trivial neutron fraction 
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in the nucleon cosmic rays above the GZK cutoff. The 
magnitude of the neutron fraction depends on the energy 
of the nucleon and the sources, and it increases with the 
increase of energy, as the time dilation effect is stronger 
for neutrons with higher energy. Detailed features of the 
neutron fraction and nucleon spectrum may also provide 
information on the sources of the nucleon cosmic rays, 
such as their distance, the nucleon spectrum of the ini- 
tial sources and their neutron/proton ratios. For example, 
assuming that the nucleons with energy above the GZK 
cutoff are from a point source with an uniform nucleon 
spectrum (i.e., with a constant value of density distribu- 
tion) at a distance far away (explicit calculation shows 
that 3A P is enough), the neutron/proton ratio should reach 
the equilibrium value of A^ ff /Ap ps 0.19 at the pion pro- 
duction threshold and 0.33 at the A production threshold, 
independent of the neutron/proton ratio of the source. 
Adopting a small cross section at threshold could reduce 
A° ff /A p to around 0.05, in agreement with explicit model 
calculation [12| . Thus, there must be a significant neutron 
fraction in the nucleon cosmic rays above the GZK cutoff, 
even if the sources only emit protons. 

However, the neutron fraction is negligibly small be- 
low the GZK cutoff, since only the neutrons within a dis- 
tance comparable to the mean free path of decay l n ~ 1 
Mpc can reach the earth. The protons with energy below 
the GZK cutoff may interact with the relic photons via 
electron-position pair production, but the energy loss is 
small, of only 0.1% compared to 20% for pion production. 
Therefore such protons may come from any source with 
distance within as large as half the Hubble length 0. We 
thus expect a small neutron/proton ratio below the GZK 
cutoff. 

For the Z-bursts of particle productions, there should 
be a symmetry between the nucleon and antinucleon pro- 
ductions in the standard model. Thus there should be 
equally numbers of protons and antiprotons in the nu- 
cleon cosmic rays from the Z-bursts. Balloon and satel- 
lite investigations can provide direct measurements of the 
cosmic rays of proton and antiprotons respectively. Un- 
fortunately, the available measurements can only reach to 
the energy scale of 10 15 eV still 4^5 orders be- 
low the energy region for our purpose. One of the meth- 
ods to identify the hadron species in the cosmic rays is 
by measuring the muon charge ratio of the air 

shower by the primary cosmic rays. It has been known 
that the / n~ ratio may provide information on the neu- 
tron/proton ratio in the primary cosmic rays |14j . With 
the improvements in knowledge of particle productions in 
hadron-hadron interactions from accelerator experiments 
and also in understanding of particle production and prop- 
agation mechanism in the atmosphere, the ratio, 
in combination with other coincidence measurements, can 
reveal the species in the primary beams. However, we will 
show that the ratio should be similar for primary 

neutron and antiproton beams. Thus the fi + / fj,~ ratio is 
hard to distinguish between the neutron and antiproton in 
the nucleon cosmic rays. With the above prediction that 
the neutron fraction is negligibly small below the GZK 



cutoff, we suggest to use the [i + / fi ratio as a measure- 
ment of the antiproton fraction in the cosmic rays below 
the GZK cutoff. 

In fact, the muons in the air showers are mainly from 
decays of pions and kaons produced in the interactions 
of the primary cosmic rays with the atmosphere |14II15| . 
The very high energy secondary pion and kaon cosmic 
rays can be considered as from the current fragmentation 
of partons in deep inelastic scattering of the primary cos- 
mic rays with the nucleon targets of the atmosphere in a 
first approximation |16) . We also consider only the favored 
fragmentation processes, i.e., the tt + , which is composed 
of valence u and d quarks, is from the fragmentation of u 
and d quarks in the nucleon beam, and the 7r _ , which is 
composed of valence u and d quarks, is from the fragmen- 
tation of u and d quarks [17} ■ Similarly, the K + ', which is 
composed of valence u and s, is from the fragmentation of 
u and s quarks, and the K~ , which is composed of valence 
u and s, is from the fragmentation of u and s quarks. The 
/i + is from the decay of a tt + or a K + and the [x~ is from 
the decay of a ir~ or a K~. We can roughly estimate the 
muon charge ratio by 

/j+ = Jo dx { [u(x) + d{x)] + k [u(x) + s(x)} } 
M~ dx {{d(x) +u(x)} + k[u(x) + s(x)}} ' 

where q(x) is the quark distribution with flavor q for the 
incident hadron beam and K^0.1^0.3isa factor re- 
flecting the relative muon flux and fragmentation behavior 
of K/ir. Secondary collisions do not influence the above 
estimation but can significantly reduce the final detected 
muon energies, since the current parton beams still keep 
their flavor content and act as the current partons after 
the strong interactions with the partons in the atmosphere 
targets. Adopting a simple model estimation of the parton 
flavor content in the nucleon without any parameter |18j. 
we find that ~ 1.7 for proton and /i + ~ 0.7 

for neutron. This simple evaluation is in agreement with 
the empirical expectation of /i + / ss 1.66 for proton and 

ps 0.695 for neutron as well as that in an ex- 
tensive Monte Carlo calculation thus it provides a 
clear picture to understand the dominant features for the 
muon charge ratio by the primary hadronic cosmic rays. 
For the ratio for antiproton, it is equivalent to 

the ratio for proton by using eq. (|17f) . thus we find 

~ 0.6 for antiproton, which is close to that for neu- 
tron. The [i + / fj,~ ratio for antineutron is also equivalent 
to the ratio for neutron, and it is ~ 1.4, 

which is close to that for proton. It is hard to distinguish 
between the primary neutrons and antiprotons (or protons 
and antineutrons) by the ratio of the air shower, 

unless very high precision measurement is performed and 
also our knowledge of the muon charge ratio for each nu- 
cleon species is well established. 

For the nucleon cosmic rays with energy above the 
pion photoproduction, there is an admixture of neutrons 
and antinucleons with the protons. It is hard to make a 
clear distinction between the neutron and antiproton (pro- 
ton and antineutron) by the muon charge ratio fi + / ji" . 
Also the number of nucleons should be very limited by 
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the GZK cutoff suppression. But there will be only pro- 
tons and antiprotons in the cosmic rays below the GZK 
cutoff from our above discussion. With the collection of a 
sizeable number of events with enough statistics, we ex- 
pect to detect the antiproton content in the cosmic rays 
by the measurement of the /i + / fi~ ratio. The difference in 
the muon charge ratios between protons and antiprotons 
will be enhanced if there is an energy cut on the detected 
muons, as the ratio \u(x) + d (xj\ / [d(x) +u(x)] (i.e., the 
dominant contribution in eq. I|17|l) has a strong x depen- 
dence to increase for proton and decrease for antiproton 
at large x. There could be also a composition of heavy nu- 
clei in the ultrahigh energy cosmic rays [23j- The nucleus 
beams are estimated, by using eq. (JTJJ, to have the muon 
charge ratio fi + / ji~ ~ 1.1, which is different from both of 
these for protons and antiprotons. However, methods have 
been developed to identify the primary nucleon and nu- 
clei by other information such as Cerenkov radiation |21| . 
The secondary muon size distributions have been also used 
in available experiments to identify the mass composition 
of primary cosmic ray extensive air showers (22123) . We 
suggest to add muon charge as a further information to 
identify the antiproton composition. In combination with 
other information, it is thus possible to detect the hadronic 
composition by measuring the muon charge ratio for the 
cosmic ray events with ultrahigh energy. 



We need to point it out here that none of the avail- 
able facilities for ultrahigh energy cosmic rays is aimed at 
measuring the muon charge ratio as their physical moti- 
vation, so it is not practical to expect the measurement of 
the muon charge ratio for the ultrahigh energy comic ray 
showers within a short period. In principle, there should be 
no difficulty to measure the muon charge ratio by the avail- 
able techniques of detectors |23], but this feature needs to 
be incorporated into design of new generation facilities for 
the physical goal we suggested. It is also not practical to 
measure the muon charge ratio on an event by event ba- 
sis. The reason is that not all particles in an ultrahigh 
energy cosmic shower can be detected: the detectors are 
distributed in arrays, and can not cover all region. So some 
muons should be lost in the detection procedures. There- 
fore the muon charge ratio should be studied from a large 
number of events in a sense of statistics. Also for our goal 
of identifying the nucleon species by their muon charge 
ratio, the muons should be detected at distances near the 
cores of the cosmic ray showers, and the energy of the 
muons should be high enough to guarantee the detected 
muons of being from the cascading fragmentation of the 
current partons in the primary cosmic rays in consequent 
multi-scattering processes. We should point out that there 
are technical limitations of the highest energy of muons 
to be detected, so we should not expect that the detected 
muons are from the the first fragmentation of the partons 
in the primary cosmic rays, but rather the current partons 
in the primary cosmic rays should still keep at the lead- 
ing current partons in the subsequent scattering processes 
so that their flavor information is not completely lost at 



the highest energy of muons to be detected . We should 
also point out that the flavor structure of the partons is 
not explicit included in the available Monte Carlo simu- 
lations, as the physics picture of the current fragmenta- 
tion we described above has not been incorporated in the 
available Monte Carlo codes or models. In similar to the 
available air shower profiles with muon size distributions 
|22| . detailed features with also muon charge and energy 
distributions for the air shower profiles of primary pro- 
ton and antiproton cosmic rays are needed and waiting 
to be constructed by further explicit researches. This can 
be only realized by new generation experiments aiming 
at detecting the antiproton content in the cosmic rays, by 
collections of a large number of events of cosmic ray show- 
ers, together with their detailed muon charge and energy 
information (perhaps with also orientation distribution in 
the earth's magnetic field) contained in the measurements. 

Our results not only work for the Z-burst hypothe- 
sis, they also apply to any model involving the neutron 
and antiproton contents in the cosmic rays. The cosmic 
rays above and the GZK cutoff can be attributed as origi- 
nated from the decays of relic metastable superheavy par- 
ticles clustered as dark matter in the galactic halo 
or from the topological defects that are left over from 
the early-universe phase transitions caused by the spon- 
taneous breaking of symmetries [Hj ■ There should be also 
some composition of neutrons and antinucleons in the cos- 
mic rays in these models. In fact, the role played by the 
time dilation effect of the neutrons has been discussed in a 
number of literature pTll25j . But the combination with the 
proton/neutron changing properties in their propagation 
makes this dilation effect more significant. We thus obtain 
an interesting scenario that the proton and neutron with 
energy above the GZK cutoff change between each other 
upon interaction with the relic photons in their propaga- 
tion, and this indicates a non-trivial neutron fraction in 
the nucleon cosmic rays above the GZK cutoff. Our pre- 
diction of the propagation of neutrons and the suggestion 
of using the muon charge ratio to detect the antiproton 
fraction of the cosmic rays below the GZK cutoff, are ap- 
plicable to all these models. 

In summary, we studied the role of neutrons in the cos- 
mic rays by taking into account the time dilation effect and 
the proton/neutron changing properties in the pion photo- 
production by the relic photons. We predict a non-trivial 
neutron fraction in the cosmic rays above the GZK cut- 
off, no matter the nucleon cosmic rays come from sources 
with distance comparable to the GZK zone or within lo- 
cal galaxies. However, there is a negligibly small neutron 
fraction in the cosmic rays with energy below the GZK 
cutoff, unless the sources are with distance comparable to 
1 Mpc. Thus, we also suggest to use the muon charge ratio 
to detect the antiproton content of the cosmic rays below 
the GZK cutoff. This may serve to provide some more 
information concerning several existing models of using 
Z-bursts, superheavy particles, and topological defects as 

1 Or the muon detector should be placed deep underground 
for the purpose to detect the muons with initially very high 
energy 
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origins for the observed events of cosmic rays with energy 
above and near the GZK cutoff. 
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